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We find tiiat a repulsive potential of graphene in the presence of a magnetic field has bound 
states that are peaked inside the barrier with tails extending over ^(A'^ + 1), where £ and N are the 
magnetic length and Landau level(LL) index. We have investigated how these bound states affect 
scaling properties of the induced density of filled Landau levels of massless Dirac fermions. For 
chiral fermions we find, in strong coupling regime, that the density inside the repulsive potential 
can be greater than the value in the absence of the potential while in the weak coupling regime we 
find negative induced density. Similar results hold also for non-chiral fermions. As one moves from 
weak to strong coupling regimes the effective coupling constant between the potential and electrons 
becomes more repulsive, and then it changes sign and becomes attractive. Different power-laws of 
induced density are found for chiral and non-chiral fermions. 
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I. INTRODUCTION 



Non-relativistic massless Dirac electrons exist in two- 
dimensional graphene layers near K and K' Brillouin 
points [Jjll. Energy dispersions form Dirac cones with 
conduction and valence bands meeting at the Dirac point. 
The wavefunction has two components: the first compo- 
nent gives the probability amplitude of finding the elec- 
tron on A carbon atoms while the second component 
gives the amplitude of finding it on B carbon atoms. In 
the presence of a magnetic field the Dirac cones split 
into Landau levels with some unique features (^-Q. A 
Landau level with zero energy that is independent of 
magnetic field develops with chiral wavefunctions. Other 
non-chiral Landau levels of conduction and valence bands 
have opposite energies but their wavefunctions are iden- 
tical except for phase factors [9l.[To|. The energies of these 
Landau levels depend non-linearly on magnetic field 



En = sgn{N)E,nV2\Nl 



(1) 



where the energy separation between these LLs is set 
by Em = -^Ul [il. Valence band LLs are labeled 
A'^ = — 1,— 2,— 3,... with decreasing energy, while conduc- 
tion band LLs are labeled N = 1, 2, 3, ... with increasing 
energy. The zero energy LL has = 0. 

These unusual graphene LLs respond rather differently 
to presence of potentials [l3l - [l6j in comparison with LLs 
of ordinary semiconductors : Confinement and decon- 
finement transitions [isj and bound states forming in- 
side an antidot [11] through complete Klein tunneling 
are predicted. Here we consider the electron density in 
the presence of a rotationally invariant and repulsive po- 
tential V{r) with strength V and range R. The electron 




FIG. 1: We have plotted schematically electron density for a 
filled Landau level. The electron density in the absence of the 
potential is represented by the dashed horizontal line. 

In graphene LLs, under certain conditions, induced density 
can be positive (a), in contrast to the usual case where it is 
negative (b). 



density of the A^th filled LL is given by 



nAr(r) 



E 



(2) 
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where eigenstates 5'5v(f) and eigenvalues En {J) are la- 
beled by LL index N and half-integer angular momen- 
tum quantum number J. In the absence of a local- 
ized potential the dimensionless density takes the value 
PnM{r) = We define the induced density as the dif- 
ference between densities with and without the potential 
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Graphene barrier has a natural energy scale, Ec = 
It is interesting to note that the ratio between this energy 
scale and the energy scale of LLs is given by the ratio 
between two length scales of the problem: 
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J' 



(4) 
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We find that the correct scahng function of the induced 
density has the form 

fAn^{r)^SMi^,^^,j). (5) 

It is different from the scaling function of ordinary LLs 
where one would expect that appears instead of 
In graphene -J^ does not contain non-perturbative effect 
of the formation of bound states in the barrier, and it 
cannot be used instead of Moreover, the variable 
is inappropriate since it becomes infinitely large at zero 
magnetic field (note E„i — at _B = 0), which leads to the 
unphysical result that the scaling function is independent 
of 

The value of the dimensionless electron density at the 
center of the potential (r = 0) gives a good indication 
of how strong the effective coupling constant between the 
repulsive potential and electrons is. We will thus call this 
dimensionless induced density at r = 0, with sign change, 
as the effective coupling constant 

«^(;^, f ) = -^'An^.(0) = -s^{Q, |). (6) 

As one moves from weak to strong coupling regimes the 
effective coupling constant becomes increasingly more re- 
pulsive, and then starts to decrease, passing through zero, 
and becomes attractive, see FigsH] and [2j The strong, 
intermediate, and weak coupling regimes correspond to 
Y ^ 1, J ~ 1, and y ^ 1, respectively. In the strong 
coupling regime a repulsive potential can effectively at- 
tract electrons, making the induced density positive. The 
physical origin of this effect is the formation of bound 
states that are peaked inside the barrier with tails ex- 
tending over 1{N + 1). We stress that these states are 
not resonant states of the repulsive potential in graphene 
at B=o (l7l - [l9| since the extent of the wavefunctions is 
finite and their energies form discrete spectra, unlike res- 
onant states. 

In the limit V/Ec ^ 1 we find the following power 
laws: 

with (Jat = 2 for the chiraljlOl = LL, but with 5n 
for non-chiral N — 1 LL. We find that go{R/i) changes 
sign near 1. A plot of gQ{R/£) is shown in Fig (21 

This paper is organized as follows. In Sec. II we present 
probability wavefunctions of boundstates in weak, inter- 
mediate, and strong coupling regimes. Using these prob- 
ability wavefunctions we investigate scaling properties of 
induced densities for chiral and non-chiral fermions in 
Sec. III. We also compute the boundary between positive 
and negative values of induced densities in the parameter 
space. Summary and discussions are given in Sec. IV. 




FIG. 2: Y-axis represents values of go{R/£) obtained by data 
collapse for N=0 and < 0.7. Note that this quantity is 

proportional to — qjv(-|-, -f). Values {R/£, V/Ec) used in the 
data collapse are shown in the inset. Note that the induced 
density changes sign near R/£ — 1. Crosses represent the re- 
sults obtained by treating V{r) in the 2nd order perturbation 
theory. These results demonstrate that perturbative methods 
cannot be applied in the intermediate and strong coupling 
regimes. 



II. MODEL CALCULATIONS OF EXACT 
WAVEFUNCTIONS 

The electron density of a filled LL is computed from 
single electron wave functions, see Eq.Q. These wave 
functions are chosen to be eigenstates of the following 
model potential 

In the strong coupling regime states with small values of 
I J| are peaked near r = and are particularly relevant 
for the density inside of the potential. 

A. Energy spectrum 

We compute eigenenergies En {J) as a function of J by 
solving Dirac equations. We choose the magnetic vector 
potential as ^ = x, 0). The two-component wave- 

functions of massless Dirac electrons obey the following 
Hamiltonian 

H ^vpS- {p+ -A) + Vif) (9) 
c 

where the Fermi velocity is vp, the Pauli spin matrices 
are a = {o-x,ay). Eigenstates ^ifir) are also eigenstates 
of angular momentum operator 

Jz = -ide + Y^ (10) 
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where is a Pauli spin matrix and 9 is the polar angle. 
Since Jz commutes with Dirac Hamiltonian of iVth filled 
LL eigenstates must have the following form 



Kir, 9) 



^ gi(./-i/2)e 



XAir) 



(11) 



FIG. 3: Single particle energy spectrum of massless Dirac 
fermions in the presence of repulsive potential and mag- 
netic field. Energies of five Landau levels are shown for 
= —2, —1, 0, 1, 2. Dashed lines represent unperturbed ener- 
gies, and those that deviate from these are energies of bound- 
states. (At B = 20T, R = lOnm, R/g = 1.74, V = O.leV). 



with half-integer values of angular quantum number J. 
For each N allowed values of J are displayed in TABLE 
1. Some eigenenergies are shown in FiglJ] Although this 
spectrum resembles the spectrum of ordinary LLs[2l| the 
wavefunctions of the eigenstates are rather different. 



TABLE I: Eigenstates are labeled by two quantum numbers: 
LL index A'^ and angular momentum quantum number J. Pos- 
sible values of A'^ and J are listed. 



N 


J 






2 


3 113 5 
C^'^' 2^ ^' 9.^ ' ' ' 


1 


113 5 

9.' 9.y 9^ 9' " ' 





13 5 

5' 9' 2' " ' 


-1 


113 5 

9' 2' 9' " ' 


-2 


3 113 5 

9^ 9' 9' 9' 9' "' 







B. Solutions in strong coupling regime 



TABLE II: How tp^ and ipn^m are related to each other. For a given J some possible values of (n, m) are listed. Empty boxes 
indicate that (n, m) do not exist. 



J ^--^ 


(n_3,m_3) 


(n_2,m_2) 


(n_i, m_i) 


(no, mo) 


(ni, mi) 


(712, m2) 




3 

9 


(-3,1) 


(-2,0) 








(2,0) 


(3,1) 


1 

9 


(-3,2) 


(-2,1) 


(-1,0) 




(1,0) 


(2,1) 


(3,2) 


1 

9 


(-3,3) 


(-2,2) 


(-1,1) 


(0,0) 


(1,1) 


(2,2) 


(3,3) 


3 
9 


(-3,4) 


(-2,3) 


(-1,2) 


(0,1) 


(1,2) 


(2,3) 


(3,4) 



It is instructive to study properties of probability wave- 
functions in the strong coupling limit, where R/l — )■ 0. 
Eigenstates (xa,X-b) determined by the pair of cou- 



pled first order differential equations 



-idxXA + i 
'idxXB - i 



i(J+i) + ix' 



XA = t{x)XB 

Xb = <^{x)xA 



, (12) 



4 



where x = r/£ is the dimensionless coordinate and e{x) 




FIG. 4: Boundstates of strong coupling limit: 

for V/Ec = 1.82, R/£ = 0.01, R = Inm, V = 1.2eV, and 

E = 0.009eV. Inset: A (B) component is dashed (dotted) 

line. 



For r ^ R the effect of the potential is negligible, and 
solutions (l^ are 



V'f (f) = ^n.mir) = Cn 



-sgn(n)i0|„|_i^,„(f) 

'/'|n|,m(^) 



(13) 



Here n and m are integers with to > 0. We define 
sgn(n) = —1,0, 1 for n < 0, n = and n > 0. The nor- 
malization constant C„ = 1 for n = and C„ = 1/V2 for 
n 7^ 0. Applying angular momentum operator to Ea. (jl3p 
and comparing with Eq. ljlll) we find that the quantum 
numbers J and {n, ni) are related to each other through 



J 



1/2. 



(14) 



For a given value of J there are infinitely many possible 
values of (n,TO), and we will order them with an index 
i. How (J, i) are related to (n, to) is given in Table II. In 
Eq. ([T3l) the wavefunctions (j)n,m{r) are the Landau level 
wavefunctions of ordinary two-dimensional systems [2^ 



(l>n,m{r) = Ai,mexp i(n-TO 



X L 



(n+m— |m-n|)/2 I 



(15) 



where An^m are normalization constants. 

To find solutions that are valid for all r, we solve the 
Dirac equations numerically using confluent hypergeo- 
metric functions [Tsj. The obtained numerical results are 
shown of FigHl Near x > R/l the value of the eigenfunc- 
tion ^'iiif) is approximately equal to ■(/)/(_/?), given by 



Ea. (IT3)) . Both A- and B-components of the eigenfunction 
vary rapidly for r < R. For J = 1/2 the A-component 
of the wavefunction xa{x) is peaked at x = while B- 
component is peaked at a; = R/i. These peak values 
of both xa and xb approach finite values in the limit 
R/£ — > 0. However, as V/Ec increases they also increase. 
The amount of jump between and R is, consistent with 
Eq.dUD, 



^AxA = -^(XA(i?)-XA(0))«e(0) 



R/e 



dxxB (16) 



with xs(0) = and e(0) = {E - V)/Em. For J = -1/2 
the properties are the opposite to those of J = 1/2: the 
B-component of the wavefunction XBix) peaked at x = 
while A-component is peaked at a; = R/£. The amount 
of jump between and R is 



iAxB « e(0) 



R/t 



dxxA- 



(17) 



In the strong coupling regime R/£ < 1 probability 
wavefunctions can be significant inside the potential, see 
FiglS] These states are peaked inside the potential range 
R and have tails extending over the length of order 
£{N + 1). Examples of such states with N = 0,1,2, 
and 3 are displayed in FiglS] As N increases the ex- 
tend of |^'jv^^^(r)P outside the potential increases ap- 
proximately as £{N + 1). However, the strength of peak 
within the range R decreases with increasing N. We 
stress that these states are not resonant states of the re- 
pulsive potential since the extend of the wavefunctions 
is finite and their energies form discrete spectra, unlike 
resonant states. States shown in FigjS] contribute to a 
positive induced charge since the probability of finding 
an electron inside the potential has increased compared 
to the probability in the absence of the potential. This 
is a non-trivial effect of the interplay between effects of 
quantization of LL and Klein tunneling (23l. [23| . 



C. Solutions in intermediate and weak coupling 
regimes 

We show how probability wavefunctions change as R/£ 
changes from weak to intermediate coupling regimes. 
In the perturbative regime i?/£ 3> 1 the exact prob- 
ability wavefunction at r = is smaller than that of 
the unperturbed probability wavefunction, as shown in 
FiglS{a)- Note that cusps in the wavefunctions are neg- 
ligible. Figinjb) displays the exact probability wavefunc- 
tions ^^'^{t)\'^ in the intermediate regime of R/£ ^ 1. 
Both A and B components of it have cusps at r = i?. 
Note that at r = the exact probability wavefunction 
is larger than that of the unperturbed probability wave- 
function. 
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FIG. 5: Strong coupling regime: the parameters are Ec — 
0.1316 eV, Em = 0.0395 eV, V = 0.26eV, V/Ec = 1.98 and 
R/i = 0.3. (a) Solid line is for iV = and E = 0.017 eV. 
Dashed hne is for iV = 1 and E = 0.061 eV. Inset: B- 
component of \'^Q^^^{r)\'^ and A- (dashed) and B- (solid) 
components of ^''^('^01^ the absence of potential, (b) 
Solid line is for iV=2 and E = 0.083 eV. Dashed line is for 
A'' = 3 and E = 0.101 eV. Inset: A- and B-components of 
|^'^^''^(r)p and |'I'3^^^^(?')p in the absence of potential. 



III. SCALING FUNCTION OF ELECTRON 
DENSITY 



We have performed an extensive numerical evaluation 
of the electron density. The dimensionless induced den- 
sity for = filled LL, i^AnQ{r/R), is plotted in FigE] 
for various values of {R/£, For R/i <1 the induced 
density inside the barrier is positive, which is in sharp 
contrast to what usually happens in a barrier. The for- 
mation of bound states inside the barrier, as discussed 
in Sec. II B, is responsible for this effect. Note the in- 
duced density oscillates as a function of r/R. We have 
tested that the total integrated density is equal to the 
total number of electrons in the LL. The induced density 




r/R 

FIG. 6: (a) Weak coupling regime: |*o ^''^(r)p for i?/^ = 2.0, 
V/E^ = 2.172, R = llnm, V = 0.13eV, and E = O.llOeV. 
A (B) component is solid (dashed) line. Inset: dashed line is 
without potential and solid is with potential, (b) Intermediate 

coupling regime: 1*0 -R/^ = 0-6, V/Ec = 0.607, 

R = 5nm, V = 0.08eV, and E = 0.013eV. 



satisfies the following two-parameter scaling function 

e^Ano{r/R) = so(r/i?, ^,R/e). (18) 

As R/£ increases the sign of Ano(O) changes sign from 
plus to minus, which is shown in FiglT] Physically this 
means positive induced density changes to negative in- 
duced density. For < 0.7 our numerical results dis- 
play data collapse, and is consistent with the following 
power law 

I'AnoiO) « {^rgo{R/£), (19) 

see Fig 12] For a larger range of < 2.5 an approximate 
data collapse can be obtained, see FigEl 

The second order perturbative calculation in V{r) 
agrees with scaling result when R/£ > 2.5, but disagrees 
when R/£ < 2, as shown in Figl^l It is noteworthy 
that go{R/£) takes the minimum value near R/£ 1.5. 
This implies that the electron density is depleted most 
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strongly for R/i « 1.5. However, further decrease in 
has the opposite effect of increasing more penetra- 
tion of electrons into the barrier. Near R/l — 1 the 
scaling function go{R/£) changes sign. For R/£ < 1 elec- 
trons accumulate in the barrier and the density becomes 
greater than the density of the unperturbed LL. This de- 
pendence on R/£ is thus strongly non-linear. For larger 
values of V/Ec the boundary between positive and neg- 
ative induced densities is displayed in Fig|9l 
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FIG. 9: Ai" = LL states must not overlap with A*' = ± 1 
LLs. Induced density at r = is positive for Rji < 1. 



and R/£ < 1. (f , -|-) = 
(1, 1.215) (solid), (f , ^) = (0.60, 0.607) (dashed), (f , -^) = 



FIG. 7: (a) rn(r) for iV 

V 

-Be 



(0.30,0.334) (dot). Dashed-dot line is e'^n{r) in the ab- 
sence of the potential. (b) For A*' = and R/£ > 1. 
if^ wj = (1-74,2.582) (solid), (f , -|-) = (1.50,1.823) 
(dashed), (f , -|-) = (1.22, 1.063) (dot). 
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FIG. 8: Approximate data collapse of go{R/£) is obtained for 
a larger range of ^ than the one used in Fig[2] . Values of 
{^,R/£) used are shown in inset. 




FIG. 10: £^n{r) for A^ = 1. {f,^) = (1.5,3.038) (dashed), 
{y,wJ = (0.6,0.607) (dot), (f ,-|:) = (0.01,0.607) (dashed 
dot). Solid line is £^n(r) in the absence of the potential. 



The induced density for iV = 1 filled LL shown in 
FigfTOl satisfies a similar scaling relation as that of = 
LL 



£''^ni{r/R) = si(T/R,^,R/£) 

V 



(20) 



At the origin r = we find, for -g- ^ 1, the following 
scaling result 



£'Anm^{^)''9i{R/e), 



(21) 
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FIG. 11: gi{R/£) is obtained by data collapse for = 1 and 
^ ^ 1 . Values oi R/£) used are shown in inset. 
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FIG. 12: N ^ 1 LL states must not overlap with iV = 0, 2 
LLs. The range of R/£, where the induced density at r = is 
positive increases with increasing V/Ec- 



where Si — I. This result is obtained in the range 
V/Ec < 0.1 by data collapsing numerical data points, 
see FigfTTl The dependence of the induced density on 
R/£ is again strongly non-linear: gi{R/£) takes the min- 
imum value near R/£ = 0.7. For larger values of the 
boundary between positive and negative induced densi- 
ties is displayed in FigHH As increases the range of 
R/£ where the induced density is positive expands. Note 
that for small values of R/£, for example 0.01, the value 
of £ at B=1T is 257 A and R becomes comparable to the 
lattice constant so that Dirac equations breakdown. In 
this case smaller values of B must be used so that R can 
take larger values. 

It can shown from perturbative analysis that, for chi- 
ral fermions, the first order correction of V{r) absent, 
but the second order correction is present and is nega- 
tive. The absence of the linear terms in -g^ for the chiral 
iV = LL is a consequence of the symmetric proper- 
ties of conduction and valence band LLs. However, for 
non-chiral fermions the first order correction is present. 
These results are consistent with non-perturbative scal- 
ing results given by Eqs. ([T9|) and ([2T|) . 

We now mention some general properties of the in- 
duced density of LLs. It has a critical point yc.i, where 

dsMO^ = (22) 

dy v=yc.i 

(x = 1^ and y = f ). We find that y^.i ^ 1- The 
scaling function takes the global minimum at yd- as y 
decreases the induced density becomes most negative at 
y = yc,i and then it increases, in contrast to the lowest 
order perturbative result in V, which suggests that it be- 
comes increasingly more negative as y decreases. Below 
this critical point perturbative methods are inapplicable, 
and it separates strong and weak coupling regimes. The 
sign of the induced density changes at the second critical 



point ?/c,2, where 

SN{0,x,y)\y=y^^=0. (23) 

As y decreases below y — yc,i the induced density 
changes sign at y = yc.2 and begins to take positive val- 
ues. 

IV. SUMMARY AND DISCUSSIONS 

We find that a repulsive potential of graphene has 
bound states that are peaked inside the barrier with tails 
extending over £{N + 1). The properties of these bound- 
states change as R/£ varies, and affect the induced den- 
sity of filled LLs inside the barrier in a non-trivial way 
as a function of R/£. As R/£ decreases the induced den- 
sity inside the barrier becomes more negative, but as it 
reaches a critical value the induced density reaches an 
extremum value. Upon further decrease of R/£ the value 
of the induced density reaches zero, and, after this, it be- 
comes positive. These changes are strongly non-linear in 
R/£, and one moves successively from weak, intermedi- 
ate, and strong coupling regimes as R/£ decreases. The 
condition V/Ec ^ 1 is not sufhcient for treating V{r) 
perturbatively, and, in addition to this, one must require 
R/£ ^ 1 for both chiral and non-chiral fermions. 

For filled LLs electron-electron interactions may be 
approximated well by a Hartree-Fock (HF) methodpsj. 
The electron density in HF method is given by the sum 
of unrenormalized single electron probability wavefunc- 
tions, just as in Eq.®. So our calculation of the in- 
duced density is actually a HF result. However, our sin- 
gle electron energies are not the renormalized HF result. 
This will somewhat affect our numerical estimate of the 
boundary between overlapping LLs of Figs 151 and [T^ The 
discontinuity in the potential of Eq.® can couple states 
in K and K' valleys, which is ignored in our approach. 
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However, our ti ght- binding calculations show that this 
coupling is small fla|. 

There is a symmetry [l^ between repulsive and attrac- 
tive potentials V{r) and —V{r) so that the induced den- 
sities of these potentials are identical. In the presence 
of a repulsive or attractive potential, both charge accu- 
mulation and depletion occur, depending on the value of 
R/i, see Figsini and \T% The appearance of a charge ac- 
cumulation near a repulsive potential, for example, could 
be explained by introducing an attractive potential via 
the transformation V{r) — )• —V{r), but this same trans- 
formation would fail to explain charge depletion since 
electrons would pile up around the transformed attrac- 
tive potential. Thus charge depletion and accumulation 
cannot be explained simultaneously in either perspective 
of repulsive or attractive potential. In addition, the criti- 
cal points R/£ — j/c,i and yc,2, where the scaling function 



takes the global minimum and where it changes sign can- 
not be explained by application V{r) —V{r). It would 
be desirable to construct an analytic theory for them. 

Properties of the bound states of the potential bar- 
rier may be observed as follows. A localized potential 
may be created by a circular gate placed on graphene 
sheet. When this gate is sufficiently close to the edge 
of the sample coupling between bound states and edge 
states may be induced, and the transmission coefficients 
of edge states may reveal properties of the bound states. 
Also these boundstates may play an important role in 
transport and magnetic properties of graphene 26]. 

S.R.E.Y. thanks Philip Kim for valuable discussions on 
various aspects of this paper. This work was supported 
by the Korea Research Foundation Grant funded by the 
Korean Government (KRF-2009-0074470). 



[1] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, 
Y. Zhang, S. V. Dubonos, I. V. Grigorieva, A. A. Firsov, 
Science, 306, 666 (2004). 

[2] A. K. Geim and A. H. MacDonald, Phys. Today, 60(8), 
35 (2007). 

[3] T. Ando, J. Phys. Soc. Jpn. 74, 777 (2005). 

[4] A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. 

Novoselov, and A. K. Geim, Rev. Mod. Phys. 81, 109 

(2009). 

[5] V. P. Gusynin, and S. G. Sharapov, Phys. Rev. Lett. 95, 
146801 (2005). 

[6] Y. Zhang, Y. W. Tan, H. L. Stormer, P. Kim, Nature, 

438, 201 (2005). 
[7] D. L. Miller, K. D. Kubista, G. M. Rutter, M. Ruan, W. 

A. de Heer, P. N. First, and J. A. Stroscio, Science, 324, 

924 (2009). 

[8] K. I. Bolotin, F. Ghahari, M. D. Shulman, H. L. Stormer, 

and P. Kim, Nature, 462, 196 (2009). 
[9] Y. Zheng and T. Ando, Phys. Rev. B 65, 245420 (2002). 
[10] C. Toke, P. E. Lammert, V. H. Crespi, and J. K. Jain, 

Phys. Rev. B 74, 235417 (2006). 
[11] M. L. Sadowski, G. Martinez, M. Potemski, C. Berger, 

and W. A. de Heer, Phys. Rev. Lett. 97, 266405 (2006). 
[12] R. S. Deacon, K. -C. Chuang, R. J. Nicholas, K. S. 

Novoselov, and A. K. Geim, Phys. Rev. B 76, 081406R 

(2007). 

[13] H. Y. Chen, V. Apalkov, and T. Chakraborty, Phys. Rev. 
Lett. 98, 186803 (2007); G. Giavaras, P.A. Maksim, and 
M. Roy, J. Phys.: Condens. Matter 21, 102201 (2009). 

[14] S. Schnez, K. Ensslin, M. Sigrist, and T. Ihn, Phys. Rev. 
B 78, 195427 (2008). 

[15] P. Recher, J. Nilsson, G. Burkard, B. Trauzettel, Phys. 



Rev. B 79, 085407 (2009). 
[16] P. S. Park, S. C. Kim, and S.-R. Eric Yang, J. Phys.: 

Condens. Matter 22, 375302 (2010). 
[17] S.-H. Dong, X.-W. Hou, and Z.-Q. Ma, Phys. Rev. A 58, 

2160, (1998). 

[18] P. G. Silvestrov and K. B. Efetov, Phys. Rev. Lett. 98, 
016802, (2007). 

[19] A. Matulis and F. M. Peeters, Phys. Rev. B 77, 115423 
(2008). 

[20] It should be noted that a locahzed potential couples chi- 
ral LL states to non-chiral LL states, and thus makes the 
A'^ = band states slightly non-chiral. Since these states 
are dominantly chiral we will call them chiral states. In 
= 1 states both A and B components are present, and 
these states are non-chiral. 

[21] S. -R. Eric Yang and A. H. MacDonald, Phys. Rev. B, 
42, 10811R (1990). 

[22] D. Yoshioka, The Quantum Hall Effect (Springer, Berlin, 
1998). 

[23] M. I. Katsnelson, K. S. Novoselov, and A. K. Geim, Na- 
ture Phys. 2, 620 (2006). 

[24] N. Stander, B. Huard, and D. Goldhaber-Gordon, Phys. 
Rev. Lett. 102, 026807 (2009). 

[25] A. H. MacDonald, S. -R. Eric Yang, and M. D. Johnson, 
Aust. J. Phys. 46, 345 (1993); S. -R. Eric Yang, A. H. 
MacDonald, and M. D. Johnson, Phys. Rev. Lett. 71, 
3194 (1993). 

[26] B.-L. Huang, M.-C. Chang, and C.-Y. Mou, Phys. Rev. 
B, 82, 155462 (2010); H. Ohldag, T. Tyliszczak, R. 
Hohne, D. Spemann, P. Esquinazi, M. Ungureanu, and 
T. Butz, Phys. Rev. Lett. 98, 187204 (2007). 



